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ABSTRACT 

Aims. A model of jet precession driven by a neutrino-cooled disc around a spinning black hole is present in order to explain the 
temporal structure and spectral evolution of gamma-ray bursts (GRBs). 

Methods. The differential rotation of the outer part of a neutrino dominated accretion disc may result in precession of the inner part 
of the disc and the central black hole, hence drives a precessed jet via neutrino annihilation around the inner part of the disc. 
Results. Both analytic and numeric results for our model are present. Our calculations show that a black hole-accretion disk system 
with black hole mass M ^ 3.66Mq, accretion rate M ^ 0.54MqS"', spin parameter a = 0.9 and viscosity parameter a = 0.01 may 
drive a precessed jet with period P = I s and luminosity L = 10" erg s"', corresponding to the scenario for long GRBs. A precessed 
jet with P = 0.1s and L = 10'" erg s ' may be powered by a system with M ^ 5.59Mq, M ^ 0.74MoS"', a = 0.1, and a = 0.01, 
possibly being responsible for the short GRBs. Both the temporal and spectral evolution in GRB pulse may explained with our model. 
Conclusions. GRB central engines likely power a precessed jet driven by a neutrino-cooled disc. The global GRB lightcurves thus 
could be modulated by the jet precession during the accretion timescale of the GRB central engine. Both the temporal and spectral 
evolution in GRB pulse may be due to an viewing effect due to the jet precession. 
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1. Introduction 

Internal shock models are extensively discussed for gamma- 
ray bursts (GRBs) (Rees & Meszaros 1992; Meszaros & Rees 
1993; Zhang & Meszaros 2004), in which an individual shock 
episode of two collision shells gives rise to a pulse, and ran- 
dom superposition of pulses results in the observed complex- 
ity of GRB light curves (e.g., Daigne & Mochkovitch 1998; 
Kobayashi et al. 1999). The observed flux rapidly increases in 
the dynamic timescale of two shell collision, then decays due to 
the delayed photons from high latitudes with respect to the line 
of sight upon the abrupt cessation of emission after the cross- 
ing timescale, shaping the observed fast-rise-exponential-decay 
(FRED) pulses. However, some well-separated GRB pulses 
show symmetric structure, and their peak energy of the vFy spec- 
trum (Ep) traces the lightcurve behavior (Liang & Kargatis 1996; 
Liang & Nishimura 2004; Lu & Liang 2009; Peng et al. 2009). 
Both the temporal and spectral properties of these symmetric 
pulse are difficult to be explained with internal shocks. In ad- 
dition, the observed isiso - Ep relation (Amati et al. 2002) or 
Liso - Ep relation (Wei & Gao 2003; Liang et al. 2004; Yonetoku 
et al. 2004) also challenge the intemal shock models (e.g., Zhang 
& Meszaros 2002). 

Quasi-periodic feature observed in some GRB light curves 
motivated ideas that the GRB jet may be precessed (Blackman 
et al. 1996; Portegies Zwart et al. 1999; Portegies Zwart & 
Totani 2001; Reynoso et al. 2006; Lei et al. 2007). It is gener- 
ally believed that the progenitors of short and long GRBs are the 
mergers of two compact objects (Eichler et al. 1989; Paczynski 
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1991; Narayan et al. 1992; see recent review by Nakar 2007) 
and core collapsars of massive stars (Woosley 1993; Paczynski 
1998; see reviews by Woosley & Bloom 2006), respectively. 
Although the progenitors of the two types of GRBs are differ- 
ent, the models for their central engines are similar, and essen- 
tially all can be simply classed as a rotating black hole with 
a rapidly hyper- accreting process of a debris torus surround- 
ing the central black hole. Such a black hole-disk system drives 
an ultra-relativistic outflow to produce both the prompt gamma- 
rays and afterglows in lower energy bands. The most popular 
one is neutrino dominated accretion flows (NDAFs), involving 
a black hole of 2 ~ lOM© and a hyper-critical rate in the range 
of 0.01 ~ IOMqS-' (Popham et al. 1999; Narayan et al. 2001; 
Kohri & Mineshige 2002; Di Matteo et al. 2002; Kohri et al. 
2005, 2007; Lee et al. 2005 ; Gu et al. 2006; Chen & Beloborodov 
2007; Liu et al. 2007, 2008, 2010; Kawanaka & Mineshige 2007; 
Janiuk et al. 2007). The different direction of angular momentum 
of two compact objects and the anisotropic fall-back mass in col- 
lapsar may conduct precession between black hole and disc. In 
this scenario, the inner part of the disc is driven by the black 
hole during the accretion process. The differential rotation be- 
tween the inner and outer parts may result in precession of the 
inner part of the disc and the central black hole, hence drive a 
precessed jet produced by neutrino annihilation around the inner 
part of the disc, forming an S- or Z-shaped jet as observed in 
many extragalactic radio sources (see, e.g. Florido et al. 1990). 
A tilted accretion disc surrounding a black hole would also make 
the precession of the black hole and result in an S-shaped jet as 
observed in SS 433 (Sarazin et al. 1980; Lu 1990; Lu & Zhou 
2005), although the angle between angular momentum of black 
hole and disc is small due to that evolution of a two compact 
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Fig. 1. Schematic picture of a precessing system. 



object system may decrease the angle between them in merg- 
ers or the anisotropic fall-back mass cannot produce large angle 
between black hole and fall-back mass in coUapsars. 

In this paper, we propose a model of jet precession driven by 
a neutrino-cooled disc around a spinning black hole in order to 
explain the temporal structure and spectral evolution of GRBs. In 
our model, the global profile of a GRB lightcurve may be mod- 
ulated by the jet procession. The temporal structure and spec- 
tral evolution may signal an on-axis/off-axis cycle of the light 
of sight (LOS) to a precessed jet axis, as proposed by some au- 
thors to explain the nature of low luminosity GRBs 980425 and 
031203 (Nakamura 1998; Eichler & Levinson 1999; Waxman 
2004; Ramirez-Ruiz et al. 2005) or to present a unified model 
for GRBs and X-ray flashes (Yamazaki et al. 2003) and the ob- 
served spectral lag in long GRBs (Norris 2002; Salmonson & 
Galama 2002). 

We present both analytic and numerical analysis for jet pre- 
cession driven by a neutrino-cooled disc around a spinning black 
hole in sections 2 and 3. Simplifying the jet emission surface as 
a point source, we demonstrate the profile and evolution of a 
GRB pulse in section 4. Conclusions and discussion are shown 
in section 5. 



2007). The precession rate Q of the central black hole and the 
irmer disc is given by Q = 2GJ{r)j(p-r' (Sarazin et al. 1980). 
Note that regions with r ^ in the disc should contribute to the 
precession. The Q decreases as r increases, so one cannot expect 
a period behavior in an observed light curve from our model. 
With the continuity equation 



M = -InrLvr, 

the precession period P then can be expressed as 

2n 



P = 



(1) 



(2) 



where M is the accretion rate and the radial inflow veloc- 
ity of the disc material. Assuming that the angular velocity is 
approximately Keplerian, we have = rCl^, and the vertical 
scale height of the flow can be written as // = Cj/Qk^ where 
= (GM/r^y^^ and c, are the Keplerian angular velocity 
and the sound speed, respectively. The Vr can be estimated as 
Vr ~ aCsiH/r) (Kato et al. 2008), where a is the constant viscos- 
ity parameter of the disc. Substituting the expressions of and 
Vr into equation (3), we have 



1.42xl0^a5a5 — -] - 



(3) 



2. Model 

An accretion disk would be warped by its precession (Sarazin 
et al. 1980). We consider a spinning black hole surrounding a 
tilted accretion disc that its rotation axis is misaligned with that 
of the black hole, as shown in Fig. 1. Its angular momentum is 
J, = aOM^/c, where M is the black hole mass and a (0 < a < 1) 
is the dimensionless specific angular momentum. Since dJ = 
Inr^'Lv^dr for a ring at radius r in the disc with width dr, we get 
J{r) - dJ/d(ln r) — InPliV^, where E and are the disk surface 
density and rotational velocity. Due to the Lense-Thirring efl'ect 
(Lense & Thirring 1918), the disc material interior to a critical 
radius r^,, which is defined as /(r^) = 7*, will be aligned with 
the equatorial plane of the black hole. The outer portion of the 
disc (r ^ rp) with sufficiently large angular momentum keeps A = 1 
its orientation. This makes the black hole precess along with the 
inner disc (Bardeen & Petterson 1975). A jet dominated by the 
ejections of neutrino annihilation around the inner part of the 
disc thus would be precessed (Popham et al. 1999; Liu et al. B - I 



It is found that P is sensitive to a, M and H/r. These parameters 
are time-dependent in the GRB phenomenon; hence the preces- 
sion period should evolve with time. Regions with r ^ r^, in the 
disc should contributes to the precession, and the evolution of 
the hyper-accretion process would make the mass and the angu- 
lar momentum of black hole increase (Belczynski et al. 2008; 
Janiuk et al. 2008), hence make evolution of the precession pe- 
riod. In addition, the nutation in the accretion system even makes 
the observed profile be much complicated. Therefore, one can- 
not expect clear period information in the GRB lightcurves. If 
the periods are shorter than the accretion timescale, the observed 
lightcurve may be composed of some pulses. Occasionally, the 
lightcurves may show quasi-periodic feature, such as that ob- 
served in BATSE trigger 1425 (Portegies Zwart et al. 1999). If 
the periods are longer than the accretion timescale, the global 
lightcurve may be a broad pulse. This is diff'erent from that for 
SS433 (Sarazin et al. 1983), in which it is assumed that the pre- 
cession periods are shorter than the viscous timescale without 
rapidly evolving with time, hence the periodic lightcurve is a 
natural consequence in SS 433. 



3. Numerical Results 

The Eq. (3) show an explicit dependence of P to a, a, M, M, 
and H/r. However, the thickness of NDAF also depends on M, 
M, a and a. Similar to P, the injected neutrino annihilation lu- 
minosity L is also a function of these parameters. In order to 
illustrate both the dependences of P and L on these parameters, 
we present numerical calculation with the method by Riffert & 
Herold (1995). This method defines general relativistic correc- 
tion factors to simulate the precession period related to the spin 
of a black hole. They are written as 



2GM laGM\^ 



3GM 



+ 2a 



/aGMxi 



(4) 



(5) 
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Fig. 2. Illustration of numerical results for P as a function of M 
using different parameter sets as marked in the plot. 



Fig. 3. M as a function of M for different parameter sets as 
marked in the plot. The viscosity parameter is adopted as a = 
0.01. The Two interaction points A nd B indicate estimates of M 
and M for a given set of (P, L, a, a). 
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—dx, 



(6) 



(7) 



where r,,, , is the inner boundary of the disc. The equation of con- 
servation of mass remains valid, while hydrostatic equilibrium in 
the vertical direction leads to a corrected expression for the half 
thickness of the disc (Riffert & Herold 1995), 



.3 a.fiU 



(8) 



where Cs = ip/p)^^^, p and p are the total pressure and density 
of the disk, respectively. The viscous shear Tnp is also corrected 
as 



(9) 



and the angular momentum equation can be simplified as (Riffert 
& Herold 1995, Lei et al. 2009) 



4nH\ 

The equation of state is 

P = Pgas + Pad +Pt+ Py, 



(10) 



(11) 



where pg.^^, p^^^, Pe, and py are the gas pressure from nucle- 
ons, radiation pressure of photons, degeneracy pressure of elec- 
trons, and radiation pressure of neutrinos, respectively (see, e.g. 
Di Matteo et al. 2002; Liu et al. 2007). The energy equation is 

written as 



Gvis — Gadv + Gphoto + Gv> 



(12) 



where 2vis> Gadv, Gphoto and Qy are the viscous heating rate, the 
advective cooling rate, the cooling rate due to photodisintegra- 
tion of a-particles and the cooling due to the neutrino radiation, 
respectively (see, e.g. Di Matteo et al. 2002; Liu et al. 2007). The 
heating rate Qyis is expressed as 



3GMM D 



(13) 



The equation system consisting of Eqs. (1), (2), (4)-(13) is 
closed for an unknown precession period P. It can be numeri- 
cally solved for a given parameter set of M, M, a, and a. We 
show P as a function of M for the parameter sets (a = 0.9, 
M = 3Mo, a = 0.01), (fl = 0.9, M = 3Mo, a = 0.1), (a = 0.1, 
M = 3Mo, a = 0.1) and (a = 0.1, M = 3Mo, a = 0.01) in Fig. 
2. It is found that P varies from tens of milliseconds to 10 ks, 
if M = 0.01 ~ 10 Mq/s, a = 0.01 ~ 0.1, and a = 0.1 ~ 0.9. 
It can approach the timescale of lightcurve or be longer than the 
accretion timescale whose provide a couple of possibilities of 
lightcurve. In coUapsar scenario, the central black hole would be 
rapidly rotates, i.e., a > 0.9. For the compact object mergers, the 
spin of the black hole is not strictly as high as that in the col- 
lapsar scenario (e.g., van Putten et al. 2001). Assuming that the 
global GRB Ughtcurves are modulated by the jet precession dur- 
ing the accretion timescale of the GRB central engine, the profile 
of a pulse duration may be comparable to P. Statistical analysis 
shows that the typical durations of long and short GRB pulses 
are ~ 1 and 0.1 second, respectively (Liang et al. 2002; Nakar 
& Piran 2002). From Fig. 2, we find that the case of (a = 0.1, 
a = 0.9) yields a P value much larger that 1 seconds in the range 
of M = 0.01 ~ IOMq. For the case of (a = 0.01, a = 0.9), we 
get f = 0.1 ~ 1 seconds for M - 0.4 ~ lOMg. This is consistent 
with the observed pulse durations for long GRBs. In order to ex- 
plain the duration of short GRB pulses, our model requires lower 
a and a as well as higher M than that for the long GRBs, indi- 
cating that the short duration GRBs may be powered by much 
violent accretion process than the long ones. 
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Fig. 4. Predicted flux F (the solid line) and £p(the dashed line) with our model for a symmetric pulse (panel a)and a FRED pulse 
(panel b) with comparisons to the observations(insets). 



The observed luminosity of prompt gamma-rays may also 
place constraints on our model parameters. We assume that 
the observed gamma-ray luminosity is comparable to the in- 
jected neutrino annihilation luminosity L. Similar to f = 
P{M, M, a, a), L is a function of M, M, a, and a, written as 
L = L{M, M, a, a). It can be calculated following the approach of 
Ruff'ert et al (1997), Popham et al (1999), Rosswog et al (2003), 
and Liu et al (2007). Since the calculation of L (or P) as a func- 
tion of these parameters is a huge time-consumed task, we per- 
form our calculation for typical L values only and present our 
results with M as a function of M, a, a for a given L (or P). We 
take L = lO^i erg s-^ for long GRBs and L = lO^" erg s'^ for 
short GRBs. Based on our analysis above, we also calculate M 
as a function of M for P = Is and P = 0. 1 s in case of parameter 
set (L, a, a)=ilO^^ erg s"', 0.9, 0.01) and (L, a, aMlO^° erg s-\ 
0.1, 0.01). We show M as a function of M for different parame- 
ter sets in Fig. 3. It is found that for a given luminosity, M as a 
function of M greatly depends on the rotation of the black hole 
(see the dotted lines in Fig. 3). The accretion rate M does not 
significantly increase with M for a = 0. 1 . However, it rapidly in- 
creases with M for a - 0.9. The behavior of the function M(M) 
varies for different P (see the solid lines in Fig. 3). As shown in 
Fig. 3, M slightly increase with M if P = 1 second. However, it 
even decreases with M for P = 0.1 second. This conflicts with 
that shown in the expUcit form of P (see Eq. 3). The reason is 
that the thickness of NDAF depends on M, M, a and a (see Fig. 
9 in Liu et al. 2007). One cannot expect an explicit dependence 
between M and M for a given P. The intersections between the 
lines are estimates of M and M for a given set of {P, L, a, a). 
We find fliat M ^ 3.66Mo and M ^ 0.54MoS-' for P = 1 s, 
L = 10^' erg s"', and a - 0.9, corresponding to the scenario 
for long GRBs. For P = 0.1s, L = 10^" erg s and a = 0.1, 
we get M ^ 5.59Mg and M ^ 0.74MoS corresponding to the 
scenario for short GRBs. These solutions are generaUy consis- 
tent with the requirements of GRB productions in simulations 
for collapsars (e.g. MacFadyen & Woosley 1999) and for binary 
coalescence of a neutron star and a black hole or two neutron 
stars (e.g. Kluzniak & Lee 1998). 



4. Temporal Profile and Spectral Evolution of a GRB 
Pulse from a Precessing Jet 

As discussed above, in the framework of our model one can- 
not expect period information from the observed lightcurves 
since the precession period is time-dependent. Since the period 
is a function of some time-dependent parameters as mentioned 
above, temporal profile and spectral evolution of pulses in GRB 
lightcurves may be a direct information of jet precession since 
the jet precession may conduct an on-axis/off-axis cycle during 
a precession period for a given observer. 

As discussed in Section 1, the E^-tracing-flux spectral evo- 
lution feature is observed in some GRB pulses (e.g. Liang et al. 
1996; Peng et al. 2009; Lu & Liang 2009). The profiles of these 
pulses are generally FRED, and occasionally are symmetrical. 
These temporal and spectral features can be explained with our 
model. We just illustrate the lightcurve and the spectral evolution 
for a point source with arbitrary radiation intensity in the axis 
with an arbitrary precession period with ultra-relativistic veloc- 
ity in the jet axis. As shown by Granot et al. (2002), assuming 
the emitting region as a point source in the jet axis, the calcula- 
tion can give reasonable results without any assumption on the 
jet structure. Therefore, we adopt the point source assumption in 
our calculations. We just illustrate the lightcurve and the spectral 
evolution for a point source with arbitrary radiation intensity in 
the axis with an arbitrary precession period for an observer (on- 
axis and ofl'-axis) at the rest frame in Section 4. If the emitting 
region is a shell of the jet with certain opening angle, the peak 
of the pulse would be flattened in case of a uniformed jet. Our 
calculation is followed by that present in Granot et al. 2002. 

The observed flux F and Ep would be amplified due to the 
Doppler effect, F = F ^(1 -j3f /(I -/Scosl'y, Bp = £^o(l-j6)/(l- 
yScos*!*). The observed time scale would be f = fo(l -y6cos*I')/(l - 
j8), where the subscript means the "on-axis" quantities, *P is 
the view angle between the jet axis and the LOS and yS = (F^ - 
I)''^/r, r is the Lorentz factor. From Fig. I, we have 

COS^P = COSffjetCOSeobs + sin6'jetSin6'obsCOS(0jet - <f>obs), (14) 

where i^jet = 0jet,o + 2nti/P (fi is the time in the rest frame of 
the central engine) and z-axis in Fig. 1 is the direction of angular 
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momentum of the outer part. We assume T = 300 in our calcula- 
tions. We compute the observed flux and peak energy of the vFy 
spectrum in a precession period corresponding to an observed 
pulse. Figure 4 demonstrates the initial "off'-axis" and initial "on- 
axis" lightcurves and corresponding Ep evolutions for an point 
source with arbitrary flux intensity and spectral hardness in the 
jet with arbitrary precession period P. The adopted parameters 
are Ooh. = 1°, 0obs = 90°, fljet = 4°, and ,^je,,o = 0° for ini- 
tial "off'-axis" observer and ^obs = 1°, lAobs - 90°, 0jet = 2° and 
0jet,o = 0° for initial "on-axis" observers. Two samples of the ob- 
servations are also shown in Fig. 4 for comparisons. It is found 
that our model can re-produce both FRED and symmetric pulses 
with -tracing-flux behavior, depending on the initial on-axis 
of ofl'-axis observations. 

5. Conclusions 

We have suggested that the differential rotation of the outer part 
of a neutrino dominated accretion disc may result in precession 
of the central black hole and the inner part of the disc, hence may 
power a precessed jet via neutrino annihilation around the inner 
part of the disc. Both analytic and numeric results are present. 
Our calculations show that for a black hole-accretion disk sys- 
tem with M ^ 3.66M0, M ^ 0.54MoS-\ a = 0.9 and a = 0.01 
may drives a precessed jet with P - 1 s and L = 10^^ erg s"', cor- 
responding to the scenario for long GRBs. A precessed jet with 
P = 0.1s and L = 10^" erg s"' may be powered by a system with 
M =^ 5.59Mo, M =^ 0.74MqS-\ a = 0.1, and or = 0.01, possibly 
being responsible for the short GRBs. These results are gener- 
ally consistent with simulations for long and short GRB produc- 
tions from collapsars and from mergers of compact stars. Both 
temporal and spectral features observed in GRB pulses may be 
explained with our model. 

The correlation between Eiso (or Liso) and Ep in the burst 
frame (Amati et al. 2002; Liang et al. 2004) are difficult to be ex- 
plained in the framework of intemal shock scenarios. Our model 
suggests an E^-tracing-flux behavior within a GRB pulse due to 
the on-axis/off-axis effect for a given observer, similar to that 
proposed by Yamazaki et al. (2004). The E^-tracing-flux behav- 
ior would give rise to the observed correlations between Ei^ 
(Liso) and Ep in the burst frame. 
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